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Abstract 

We explain how masses and matrix elements can be computed in lattice QCD 
using Schrodinger functional boundary conditions. Numerical results in the quenched 
approximation demonstrate that good precision can be achieved. For a statistical sample 
of the same size, our hadron masses have a precision similar to what is achieved with 
standard methods, but for the computation of matrix elements such as the pseudoscalar 
decay constant the Schrodinger functional technique turns out to be much more efficient 
than the known alternatives. 
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1 Introduction 



In recent years, lattice QCD calculations in the quenched approximation have reached a 
new quality |l|-j3|. The renormalization of many local composite operators can be treated 
non-perturbatively, and the leading discretization errors have been removed. Conse- 
quently one would now like to perform continuum extrapolations for various hadron 
masses and matrix elements in the improved theory and compare to recent results 
obtained with more standard techniques Q. In particular, the problem of a determi- 
nation of the light quark masses can be addressed with confidence, since the complete 



renormalization is known non-perturbatively [10|. 

However, continuum extrapolations require masses and matrix elements for finite 
values of the lattice spacing which are sufficiently accurate 0. Efficient methods are 
needed to obtain such precise results. Standard methods use correlation functions in 
a sufficiently large box with periodic boundary conditions fll|| . Usually one seeks to 
enhance the dominance of the low lying hadrons by using tuned hadron wave functions 
\ T% 21 1, possibly combined with variational techniques [17, 22, 23|. 



In this paper we investigate an alternative to these standard methods. Again we 
choose a sufficiently large box but impose Dirichlet boundary conditions in time, as they 
are used to formulate the QCD Schrodinger functional . We shall demonstrate that cor- 
relation functions in the Schrodinger functional are dominated by hadron intermediate 
states at large Euclidean time. Moreover, it is shown that a time extent of 3 fm for the 
box is sufficient to extract masses and matrix elements. An advantage compared to the 
standard methods is that the pre-asymptotic decay of Schrodinger functional correlation 
functions is very slow, which means that a large signal remains at large separations. 

In the following, we briefly discuss the foundation of the method (Section ^) and 
test its applicability in practice in the quenched approximation (Section |3|). We also 
attempt to quantify the efficiency compared to more standard calculations (Section 
Finally we discuss open questions as well as possible further improvements. 

2 Correlation functions at large time separations 

We now derive explicit expressions for the representation of Schrodinger functional cor- 
relation functions in terms of intermediate physical states. Throughout this section we 
assume that the lattice Schrodinger functional is defined using the standard Wilson 



action as in ref. [24|. In this situation the relations presented here hold exactly. If one 
considers the 0(a) improved theory, as we will do later, we cannot derive the equations 
given in this section directly from the transfer matrix. However, universality implies 
that the renormalized correlation functions of the improved theory and the unimproved 
theory agree in the continuum limit. Since the correlation functions considered below 
are renormalized multiplicatively, their time dependence (bare or renormalized) is given 
by the expressions derived for the Wilson theory up to lattice spacing effects. For the 
0(a) improved theory this means that all relations derived in this section are valid for 
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physical distances large compared to the lattice spacing (up to corrections of order a 2 ). 



The correlation functions considered here have been introduced before [25-27]. In 
those references the emphasis was largely on the perturbative regime, which means 
choosing small extensions of the space-time volume. By contrast, in this work we are 
interested in the correlation functions for intermediate to large volumes, i.e. extensions 
which are significantly larger than typical QCD scales. Provided that the pion mass is 
not too small, such typical QCD scales are of order 1 fm. 

The QCD Schrodinger functional is defined as the QCD partition function in a 
cylindrical geometry, i.e. periodic boundary conditions with periodicity length L in 
three of the four Euclidean dimensions, and Dirichlet boundary conditions in time at 
the hypersurfaces xq = and xq = T. Its quantum mechanical interpretation has first 
been discussed for the pure gauge theory p8| and subsequently for the theory with 



quarks [24|. In these references it was shown that the Schrodinger functional partition 



function can be written as 

Z = (f|e- TH P|i), (2.1) 

with states |i) and |f) which are given in terms of the boundary values specified at xq = 
and xo = T, respectively. In the above equation M denotes the Hamilton operator of 
QCD formulated on a torus of volume L 3 . More precisely, in the lattice theory it is 
proportional to the (negative) logarithm of the transfer matrix. 2 

Of course, the same operator describes the correlation functions when the Dirichlet 
boundary conditions are replaced by periodic boundary conditions in time. The pro- 



jector P projects onto the gauge invariant subspace of the Hilbert space [24, 28|; only 
gauge invariant intermediate states are physical and can contribute. 

For our present investigation, we have considered only the case of homogeneous 
boundary conditions, where the spatial components of the gauge potentials are set to 
zero at the boundaries and also the fermion boundary fields are taken to vanish. In this 
case we have |i) = |f) = |io) and this state carries the quantum numbers of the vacuum. 
Other choices for the boundary conditions may be of interest as well and can be treated 
similarly. 

As an example we will discuss two specific correlations, which allow for a calculation 
of the pion mass and decay constant. The generalization to other channels and other 
matrix elements is straightforward. 

We start from the dimensionless fields 

° = $ECu(y)75Cd(z), o' = pEci(yhc;w (2.2) 

y.z y.z 

and a local composite (gauge invariant) field X(x) (which will have mass dimension 
three in our applications) to define the gauge invariant correlation functions 

fx{xo) = -yW*)0>, (2.3) 



For the unimproved Wilson action, H is known to be hermitian 
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fi = -\{o'o), 



(2.4) 



where the average denotes the usual path integral average and u, d are flavour indices. 
The "boundary quark fields", Q, . . . , have been discussed in |p6| . In the lattice theory, 
Q is given explicitly in terms of the gauge fields connecting hypersurfaces xq = and 
xq = a and the quark fields on the hypersurface xq = a. Analogous properties hold 
for the other boundary quark fields. Our choices for X are X = Aq (defining, through 
eq. fl2.3|), the correlation function /a) and X = P (which gives /p), where 

A Q (x) = V , d( a; )7o75^u(a;), (2-5) 
P(x) = lp d (x)^ u (x). (2.6) 

The correlation functions fx have the quantum mechanical representation 

/x(x ) = 2- 1 ^(io|e-( T - a;o ) H PXe- :coH P|i 7r ), a < x < T - a , (2.7) 

where X is the corresponding operator in the Schrodinger picture, and the state li^) 
has the quantum numbers of the 7r + with momentum zero. To conclude that eq. (|2.7|) 
holds, one only requires that the combination of fields O, eq. has support for 

xq < a and that it carries the quantum numbers of a 7r + . The former is guaranteed by 
the very construction of the boundary fields £)C p6[| . Furthermore we have 

f 1 = 2- 1 ^(i v \e- Ta ¥\i K ). (2.8) 

It is now apparent that (for large separations xo and T — xq) the mass of the pion and 
its decay constant can be extracted. To see this explicitly we insert a complete set of 
eigenstates of the Hamiltonian, 

\n,q), n = 0,l,... , (2.9) 
M\n,q) =E^\n,q), (2.10) 

with normalization (n 1 ,q'\n,q) = <5 n ,n'^,q'- Here the energy levels in the sector of the 
Hilbert space with internal quantum numbers q are enumerated by n. Only quantum 
numbers q = ir, a shorthand for (J, P,C, I, I3) = (0, — ,+,1,1), and q = 0, which 
denotes vacuum quantum numbers are considered in the following. We do not indicate 
the momentum of the states \n,q), since both |io) and |i„-) are invariant under spatial 
translations and only states |re, q) with vanishing (spatial) momentum contribute. 

The representations given so far hold for arbitrary L, T and xq. In this paper, we 
shall be interested in the special case of the asymptotic behavior of /a(xq), f\ for large 
values of both xq and T — xq, while L remains unspecified at this stage. We include 
the first non-leading corrections but neglect any contributions which are suppressed by 
terms of order 

exp(-r J Ei 0) ), exp(-xoE[ n) — (T — xq)E { ® ) ), exp(-(T - x )£? } ) and ex.p(-xoE^ n) ) , 
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compared to the leading terms in the correlation functions. In this approximation we 
obtain 



fx(xo) ~ f p(0,0|X|0,7r)e-^ x {i + ^e"^ + ^ e - (T " xo)mG } , (2.11) 
/! « Ve- Tm -, (2.12) 



2 

where we have introduced the ratios 

(0,0|X|l,7r)(l,7rM 

lx (0,0|X|0,7r)(0 )9 r|W> ' { ' ' 

o (i |l,0)(l,0|X|0,7r) 

/x (i |0,0)(0,0|X|0,vr) ' 1 ; 

The energy difference diq = E^ — E^ is the mass of the ++ glueball and A = 
E\ — Eq is an abbreviation for the gap in the pion channel. As indicated above, we 
have dropped contributions of higher excited states which decay even faster as xq and 
T — xq become large. 

Considering the special case of /a, we find that it is proportional to the matrix 
element (0, 0| Ao|0, ir), which is related to the pion decay constant F n through 

Z A (0,0| A |0,7r) = F 7r m 7r (2m 7r L 3 )- 1 / 2 . (2.16) 

Here, Zx is the renormalization constant of the isovector axial current, and the factor 
(2m 7r L 3 ) _1//2 takes account of the conventional normalization of one-particle states (in 
our convention the experimental value of the pion decay constant is 132 MeV). 



Eq. ( 2.11 ) is used to determine m w , while the pion decay constant, F n , may be 
conveniently extracted from the ratio 

Z A fA(xo)/VTi ~ iF^m^L 3 ) 1 /^-^-^ 2 )^ 

x {l + r]le- xoA + r)° A e-( T - x °) mG } . (2.17) 

The above formulas show explicitly how masses and matrix elements can be obtained 
from Schrodinger functional correlation functions. Before we describe the numerical 
tests of the practicability, we wish to point out some properties of the present method 
and list the differences to conventional approaches. 



Equations (pTTpl^) can be expected to be rather accurate when all time separa- 



tions are larger than typical hadronic length scales (say, 1 fm, provided is not 
too small) . Schrodinger functional correlation functions decay slowly for small xq , 
leaving a large and precise signal at separations of 1 — 2 fm. This is easily seen by 
applying asymptotic freedom and a simple dimensional analysis. 
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As a comparison, consider correlation functions of standard local composite fields 
such as E x ( $t (aO$(0)>, with $(x) = ^ u (x)Ti/) d (x), &(x) = ^ d (x)j Tjoip n (x). 
At distances xq ~ 1 fm, such correlation functions are typically very small which 
usually means low statistical precision. The reason why they are small is because 
they decay like (xo)~ 3 for short time separations, as may be inferred by the same 
arguments used above. This qualitative difference arises from the fact that in 
the Schrodinger functional a dimensionless (non-local) field, J dy Cu(y)7sCd(0), is 
used to create hadronic states at the boundary. 



The ratio p, eq. ( 2.13 ) is divergent, since the state involves the bare boundary 



quark fields. However, in the final quantities of interest p is cancelled explicitly 



(see eq. (2.17)) 



• In particular, the combination Za /a(^o) /\/7i has a continuum limit for all values 
of xq . One may therefore choose some (not too large) value of the lattice spacing 
to determine the time separations xq and T — x$ where the contamination due 
to excited states is small. The same separations (in physical units) can then be 
used for other values of the lattice spacing. This holds also for the "local masses" 
which are commonly used to extract hadron masses. When one applies smearing 
or fuzzing the validity of such a statement is not immediately evident. 

• Spatial translation invariance is used fully and reduces statistical errors. 

• The present approach is similar to using "wall sources" to create hadronic states 
[jDJ, but here we keep gauge invariance in all stages of the formulation! 

3 Extraction of masses and decay constants 

In this section we demonstrate the practicability of the method in the case of quenched 
lattice QCD. 

3.1 Computational details 

We work in 0(a) improved QCD as detailed in [26|, using the non-perturbative estimates 
of the improvement coefficients c sw and ca reported in ref. [30|. 

The full 0(a) improvement of the Schrodinger functional correlation functions re- 
quires also 0(a) counterterms (with coefficients Ct,ct) at the boundary [26]. These do 
not affect the lattice spacing dependence of hadron masses and matrix elements, and 
therefore their coefficients are not very important in the present context. The only 
place where they do play a role is the size of the correction terms in eq. ( 2.11| ), since the 



lattice spacing errors in the amplitude ratios rjx^Vx are 0(a 2 ) when ct,ct are chosen 
appropriately. Otherwise cutoff effects of order a remain in these excited state contri- 



butions. We used one-loop estimates for c t ,c t [31, 3^]. For these values the 0(a) effects 
are expected to be quite small [1Q,B3]. 
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In the pseudoscalar channel we consider the correlation functions fp(xo) and f\ (xq), 
where the latter is obtained from eq. fl2.3| ) by inserting the time component of the im- 
proved axial current |25|,|(|, viz. 

(Ai^x) = A^x) + ac A \ {d; + dp) P(x). (3.1) 

The correlation function in the vector channel is defined by 

4(x ) = -a 6 ^^(m) fc (x)Cu(y)7 fc Cd(z)), (3.2) 



y 5 z 



where 



(Vi) M (a?) = ^ d (x)^ u (x) + acv%(d 1/ + d*)[xl) d (x)(T llI ,'iJj vl (x)] (3.3) 
is the 0(a) improved vector current. For the improvement coefficient cy we have used 



the values from its non-perturbative determination reported in [34|. 

We have chosen the parameters of our simulations such that L ~ 1.5 fm and T ~ 
3fm. Our calculations have been performed for four different values of the lattice 
spacing, but for the purpose of demonstrating the practicability of the method we 
restrict ourselves to results obtained at /3 = 6.0 and (3 = 6.2. These couplings correspond 
to lattice spacings a = 0.093 fm and a = 0.068 fm, when ro = 0.5 fm is used to set the 



scale [35, 36 j. For these parameters a direct comparison with results obtained using 
conventional methods @,[5],|||,|1] can De made. 

The numerical computations of Schrodinger functional correlation functions have 
been explained earlier |3(|, and further details can be found in that reference. Here we 
only mention that in addition to the previously used even/odd preconditioned version 



of the BiCGStab solver, we have also implemented SSOR-preconditioning [37,^]. The 
latter reduced the number of BiCGStab iterations needed to solve the Dirac equation 
by more than a factor 2. This turned into a gain in CPU-time of a factor of around 1.5 
in our implementation on the APE-100 machines. 

Our statistical samples consist of 1000 "measurements" of Schrodinger functional 
correlation functions at /3 = 6.0 and 800 measurements at j3 = 6.2. All statistical errors 
were computed using the jackknife method. 



3.2 Plateaux 

Let us first get a rough impression of the range of xq, where the leading term in eq. ( p. 11 ) 
dominates. To this end we follow the tradition of looking for a plateau in the effective 
mass, 

m eS (x + f ) = ^ In (f(x )/f(x + a)) . (3.4) 

Here, / denotes any of the two-point correlation functions defined above. In the pseu- 
doscalar channel we use f = f\, and the resulting effective masses are shown in Fig. |l|. 
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Figure 1: The effective mass of the correlation function f\, as defined in eq. (3.4). The three sets of 
data points correspond to (/?, k) = (6.2,0.1349), (6.0,0.1342), (6.2,0.1352) from top to bottom. In all 
cases the time extent is T « 3fm. 



There is good evidence for plateaux starting at xo ~ 1 fm and extending to approx- 
imately T — xq ~ 1 fm. As expected, the location of the plateaux is approximately 
independent of the lattice spacing. 

Turning to the vector channel, we use eq. (34) with f = kh. In this channel the 
effective mass, shown in Fig. ^, turns into a plateau only at xq ~ 1.5 fm. Furthermore, 
statistical errors grow more rapidly as xo becomes large. By comparing Fig. [I] and 
Fig. ^ one may anticipate that it is more difficult to obtain a reliable determination of 
the vector meson mass than it is to compute m n . However, this is not much different 
when standard methods are employed. 



3.3 Excited state corrections — averaging intervals 

Let us now discuss the extraction of masses and matrix elements in some detail. A 



standard method is to perform fits to the leading term in eq. ( 2.11 ). Alternatively one 
can average the effective mass over the plateau region. Both of these procedures reduce 
the error in the mass compared to just taking one point of the effective mass in the 
plateau. One must then ensure that the final statistical error is not accompanied by 
a noticeable systematic error due to a small contamination by excited states. In other 
words, one needs a rough estimate of the size of the contribution of excited states in 
the region of xq considered. 
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Figure 2: The effective mass of the correlation function k\j for a ~ 0.07 fm. 

In order to arrive at such an estimate, one requires some information about the 
values of uiq and the gap A in eq. ( 2.11]) and the analogous equation for the vector 
channel. Glueball masses from the literature [pSQ, combined with ro/a pfifl, give 



m G r « 4.3 . 



(3.5) 



In addition, we have obtained a rough estimate for A in the pseudoscalar channel by 
analyzing correlation functions with periodic boundary conditions in time, which were 
made available to us by the Tor Vergata group || . Since this analysis is not of immediate 
interest for our discussion of Schrodinger functional correlation functions we relegate the 
details to the appendix. 

The typical uncertainties associated with the gaps and the scalar glueball mass are 
of order 10% for the range of quark masses considered here. The gaps and tuq can 
now be used to take a closer look at the effective masses and obtain estimates for the 
excited state contributions and the range of xq where these are small compared to the 
statistical errors. We discuss this separately for the two channels. 



3.4 The pseudoscalar channel 

The analysis described in the appendix yields 

r A « 3.2 , 



(3.6) 



8 



0.24 



0.22 



0.2 



0.18 



i 1 1 r 



"i — i — i — i — I — i — i — i — i — I — i — r 



_i i i i_ 



_i i i i_ 



0.005 

exp(-x A) 



0.001 0.002 

exp(-(T-x ) m G ) 



Figure 3: The effective mass of the correlation function f\ against the two leading corrections. Small 
values of xo are omitted in the right half of the graph and large values of xo in the left half. Data are 
for a lattice spacing of a ~ 0.07 fm; the dashed lines indicate the agreement with linear behaviour. 



which agrees with an estimate of the same quantity using data by UKQCD Using 
eq. ( |2.11 ), the time dependence of the effective mass is given by 



m cfi (xo 



m 



m 



fj, 2sinh(aA/2) 7t -iqA _ 2 sinh(am G /2) - (T-x )m G \ 
T \ am-K 'A am, 'A J 



ma -(T-x )m G \ 
m v r /A e ] > 



(3.7) 



where the second line is valid close to the continuum limit (a A <C 1, aniQ < 1). In 
order to check whether this time dependence is reproduced by the data, we plot the 
effective masses directly against the expected form of the asymptotic correction terms, 
e~ XoA and e -(T-x )m G in Fig> | 

One observes that the data fall approximately onto straight lines. It is important to 
bear in mind that tog, A on the one hand and m e g(xo) on the other, have been obtained 
independently from correlation functions computed for different boundary conditions 



channel e t min /r (T-t max )/r 



f\ 0.2% 2.6 2.3 
/I//p, fp 0.1% 2.8 2.5 
kh 0.2% 3.0 2.2 



Table 1: Ranges of xq, where (relative) excited state contribution are smaller than e. 
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p 


K 


am n 


ref. 


ref. § 


ref. § 


6.0 


0.1338 
0.1342 


0.3529(11) 
0.3001(12) 


0.2988(17) 






6.2 


0.1349 
0.1352 


0.2430(6) 
0.2004(6) 


0.2444(9) 
0.2016(11) 


0.2007(40) 


0.2440(21) 
0.2007(26) 




P 


K 


£jbarc 




ref. § 


ref. @ 


6.0 


0.1338 
0.1342 


0.2125(7) 
0.1794(8) 


0.0943(4) 
0.0905(4) 


0.0907(8) 




6.2 


0.1349 
0.1352 


0.2165(5) 
0.1769(5) 


0.0721(3) 
0.0687(3) 


0.0727(9) 
0.0690(30) 


0.0740(35) 
0.0706(46) 



Table 2: Results in the pseudoscalar channel compared to values from the literature 



where excited state corrections have different amplitudes. The apparent compatibility 
of the data with the expected form of the correction terms is therefore non-trivial, 
and we conclude that we have a semi-quantitative understanding of the excited state 
corrections. 

It is now easy to deduce values for t m \ n and T — i max such that for t m \ n < xq < 
tm&x these corrections are below a certain margin, e, which is allowed as a systematic 
uncertainty. We list t m j n and T — t max together with the chosen values for e in Table 
Reliable estimates for are then obtained by averaging m e fl for t m \ a < xq < t ma _ x , and 
representative results are collected in Table 

A widely used method to extract hadron masses is to fit the correlation functions, 
using the expected asymptotic behaviour as an ansatz. In order to check the consistency 
of results obtained by averaging m e fr for £ m ; n < xq < t max , we have also performed single 
exponential fits to f\(xo) over the same interval. A comparison of the two methods 
shows that the estimates for m n are entirely consistent. Furthermore, the quality of 
the fits is very good, with typical values of the correlated x 2 / n df i n the range 0.6 — 0.9. 
Moreover, the stability of the fits under variations of the fitting interval has been used as 
an additional check that our values of i m i n , tmax were chosen appropriately. In Table § 
we compare our results for am^ to those obtained using conventional techniques and 
find good agreement. 

Next, we discuss the bare pion decay constant, F^ are , which is obtained from an 
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Figure 4: Influence of excited states on the ratio f\{xo) / fvixo) for a ~ 0.07 fm (open circles) and 
a w 0.09 fm (filled circles). 



average of 



i^ are = 2{m« L 3 )-V2 e (xo-T/2)m 7 



(3.8) 



over the same interval of xq, with taken from the previous analysis.^] Values for 
F^ are are also included in Table [| Note that in the improved theory which we consider, 



the renormalized decay constant is given by F n = Z^{\ + b^am^F^ |^ 

A further example for the determination of matrix elements is the combination 
m 7r F^ are /G^ are which is related to the ratio f\/ fp via 



fl 



aK x o) 



bare 



{i + vl 



e x A _|_ ^O e -(T-io)mG J 



fp(x ) G^ c {l + r# 



e -xoA + ^Og- 



(T-xo)m G 



(3.9) 



Again one may average the l.h.s. over a range of xq where excited state corrections are 
negligible. To find the proper range for the ratio f\{xo)/ fp{xo), we literally repeat the 
analysis performed before for m e ff . As shown in Fig. || the corrections of order e~ x ° A 
are of the same order as before. They originate predominantly from the denominator, 
since the PCAC relation predicts 



Vp 



A + 



~Va ■ 



(3.10) 



The leading correction terms to eq. (3.8) are suppressed by factors m^/A and m^/mG compared to 



eq. (3.7). However, we did not enlarge the interval in this case, since m n is not very small in comparison 
to the other masses. 
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Figure 5: The effective mass of the correlation function fcy against the leading correction. The meaning 
of the symbols is as in Fig. ^| 

This enhancement factor compensates for the missing factor compared to eq. ( |3.7| ), 
and a similar value of t m i n has to be chosen (see Table |l|). Results for m 7T F^ are /G^ arc 
are included in Table 0. 



3.5 The vector channel 

In analogy to the pseudoscalar case, the analysis of the correlation function ky requires 
information about the gap in the vector channel. The effective masses are plotted in 
Fig. ||. Here, the estimate for A has been obtained by tuning its value until a roughly 
linear behaviour was observed. Thus, unlike the case of the pseudoscalar, the gap is not 
known independently through a separately determined correlation function. However, 
from additional runs performed on a larger volume, we know that the contribution from 
excited states has approximately the same magnitude for L = 2.2 fm, which adds further 
credibility to the analysis of the gaps presented here. 

Our statistical errors are too large to observe a significant signal for the glueball 
contribution at large xq. The value of i max was therefore determined by requiring that 
the maximally allowed glueball amplitude be contained within the statistical errors in 
the range of x$ corresponding to exp(— (T — xo) m c) < 0.002. As before, our results for 
m p obtained using the averaging procedure were consistent with single exponential fits 
to fcy. All parameters and numerical results are listed in the tables. 

The comparison of our estimates for m p to results employing standard methods in 
Table || shows that our numbers are slightly lower, although the difference is mostly not 
statistically significant. In view of the many checks of our analysis, we are confident 
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p 


K 


am p 


ref. 


ref. § 


ref. § 


6.0 


0.1338 
0.1342 


0.508(3) 
0.480(4) 


0.487(3) 






6.2 


0.1349 
0.1352 


0.359(3) 
0.335(4) 


0.363(4) 
0.343(5) 


0.353(15) 


0.363(8) 
0.335(12) 






Table 3: 


Vector meson 


masses 





that the vector mass has been extracted correctly. It is well known that in general the 
determination of the mass in the vector channel is not easy, so that differences at the 
level of around 1.5 standard deviations are not too surprising. 

4 Numerical efficiency 

We can now assess the numerical efficiency of our method in relation to results obtained 
using conventional techniques. Comparing the errors in Tables § and ||, one has to take 
into account that the statistics for the simulation in Q is approximately the same as 
ours, whereas in refs. [|| and || the number of "measurements" is smaller by roughly a 
factor 8. If one compensates for trivial statistical factors, the tables demonstrate that 
correlation functions computed in the Schrodinger functional allow for the determination 
of hadron masses with similar precision compared to conventional methods. This is also 
the case for ratios of correlation functions like f\(xo)/ fp(xo), which serves to extract 
the combination m^F^ rc /G^ re g,|5]]. 

Another relevant issue for the overall precision is the tolerated maximum contam- 
ination by excited states, e. In order to avoid the total error to be noticeably affected 
by systematic effects, the averaging or fitting intervals must be chosen such that the 
statistical error is still significantly larger than e. It then turns out that in our approach 
one can use very small values for e without compromising the statistical accuracy. This 
is illustrated by a direct comparison to the results of ref. (8j in the pseudoscalar channel. 
From the formulas in the appendix we have 

for the analysis of ||]. Inserting our estimates for A and xp obtained from fits described 
in the appendix and the value of i m i n used in M, one obtains e 0.6% in the pseudoscalar 
channel for ref. ||. Thus, in our simulation both the statistical error and the residual 
contamination by excited states is smaller by about a factor three. 

The overall errors of the observables discussed above are similar to the ones achiev- 
able with standard methods, but perhaps - as we just argued - the Schrodinger func- 
tional correlation functions give somewhat more precise results. In addition the Schrodinger 
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functional enables one to compute the pseudoscalar decay constant with much better 
precision compared to what is usually achieved with conventional correlation functions. 
The reason for this is not entirely clear, but it may be because F^ are defined in eq. ( |3.8|) 
is obtained from a straight ratio of correlation functions times a function involving m n 
only. 

5 Discussion 

In this paper we have shown how correlation functions with Schrodinger functional 
boundary conditions can be used to compute hadronic quantities like meson masses and 
matrix elements with high precision. An integral part of our analysis was the detailed 
investigation of the influence of excited states: in the pseudoscalar channel we have used 
independent information about the gap A and the lightest glueball mass to select the 
appropriate averaging intervals. 

As explained in Section |2[ correlation functions with Schrodinger functional bound- 
ary conditions decay slowly, resulting in accurate results for masses and matrix elements. 
In particular, we have seen that our method produces very precise results for the pion 
decay constant. One may expect that a similarly good efficiency applies to pseudoscalar- 
to-pseudoscalar matrix elements such as Bk- Such applications should be investigated 
in the future. Furthermore, completely different channels like the nucleon are accessible 
and should be tried. 

All our detailed investigations have been done for L ~ 1.5 fm. How does the size 
of excited state corrections depend on L? For smaller L, we expect the dominance 
by the ground state to be similar or even better. For significantly larger L, however, 
our correlation functions might receive bigger contributions from excited states and the 
efficiency might deteriorate. We have investigated also L ~ 2.2 fm (keeping T fixed) 
for two different pairs of {(3,k). The magnitude of excited state corrections is hardly 
different from our results on the smaller system. So the location of the window of xq, 
which allows for an extraction of physical masses and matrix elements, is independent 
of L between 1.5 fm < L < 2.2 fm. Even larger values of L can not be reached with our 
computing resources but are also not necessary for the quantities studied here. 

Further improvements of the method are possible. So far we have used only the 
simplest implementation of composite boundary fields in the calculation of fi, f\ and 
ky. More refined choices of sources involving the boundary fields, such as tuned hadron 
wavefunctions, can surely be made, whilst preserving gauge invariance at all stages of 
the calculation. This might further enhance the efficiency of the method. 
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A Determination of A 

Here we describe the determination of the gap A in the pseudoscalar channel using 
correlation functions with periodic boundary conditions in all four space-time directions. 
The quark propagators were made available by the Tor Vergata group, and more details 
about the simulation can be found in || . Here we only state that A has been determined 
at P = 6.2 on a lattice of size 24 3 x 48, with the 0(a) improved action. In order to 
distinguish the correlation function computed using periodic boundary conditions from 
those defined within the Schrodinger functional , we use the letter 'C, defining 

Cpp(xo) = a 3 ^(P(x)pt( )), (A.l) 

X 

C AP (x ) = a 3 ^((^)o(x)pt( )), (A.2) 

X 

C AA (x ) = a 3 ^((^i)o(x)(4) (0)>, (A.3) 



with P(x) and (Ai)q(x) as given in eqs. (|2.6| ) and (|3.1| ). Furthermore, we here use 
P^(x) = —ip u (x)'j5\(jd{x), and similarly for A\. The spectral decomposition of the 
correlation functions in terms of the two lowest intermediate states is 

Cxx(xo) « de" IO ""x{l + xle- I » A },I = A,P (A.4) 
Cap(xo) « ^pe- lom " x {i + XAXpe -^ A }, (A.5) 

when terms proportional to exp(— m vr (T — xq)) and exp(— niG(T — xq)) can be neglected. 
We now consider the ratio 



Rap{x ) = 7^ 1—^2 ' A - 6 

[CaA x o)\ 



which - in the same approximation - is given by 

i? A p(zo) - 1 « (XP - XA) 2 e- A *°. (A.7) 
As in eq. ( |3.10 ), xp an d XA are related by the PCAC relation, 

Xp = XA, (A.8) 

mi? 



and hence 

i? A p(x ) - 1 « xl ( — ) e~ Ax °. (A.9) 



A x 2 
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Figure 6: -Rap — 1 for k = 0.1345, f3 = 6.2. The line represents a fit in the window 3a < xq < 8a. 

By fitting (R\p(xo) — 1) to the above functional form in the appropriate range of 
xq, one can extract the gap A. A typical fit is shown in Fig. || from which we obtain 
the result quoted in the text, 

r A«3.2. (A.10) 

It turns out that A depends very little on the bare quark mass, so that this result is 
used in the analysis of correlation functions in Section ^ at all values of the quark mass. 
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